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Abstract. The low energy scattering of heavy positively charged partides 
on hydrogen àtoms (H) are investigated by solving the Faddeev equations in 
configuration space. A resonant value of the pH scattering length, o = 750 ± 5 
a.u., in the pp antisymmetric state was found. This large value indicates the 
existence of a first excited state with a binding energy B=1.14xl0~ 9 a.u. 
below the H ground state. Several resonances for non zero angular momenta 
states are predicted. 



1 Introduction 

The scattering of heavy charged partides (X + ) on àtoms at kinetic energies 
smaller than the inelastic thresholds is dramatically influenced by the presence 
of atòmic electrons. Their virtual excitations in presence of the incoming charged 
particle result into long range attractive forces which dominate the low energy 
scattering. These states were proposed as a possible source of metastability in 
the pHe system Q and were found to play a determinant role in the low energy 
pH annihilation |2j . A full solution of these problems is however made extremely 
difncult by the presence of annihilation channels and only approximate solutions 
were achieved. We present in what follows a rigorous solutions for the simplest 
problem of scattering on hydrogen àtoms. 

A simple two-body approximation is first discussed in Section ||. Section [3| is 
devoted to the solution of the 3-body (X + e~p + ) problem, obtained by solving the 
Faddeev equations in configuration space. We treat with special care the pH case, 
for it exhibits the more interesting properties and constitutes moreover a realistic 
experimental challenge. Some final remarks about the present calculations and 
plans for future work are given in the conclusion. We use all along the paper 
electrònic atòmic units (m e = e 2 = h = 1). 
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2 Two-body approach 

To get a first qualitative insight into the underlying physics, it is interesting to 
consider a simple 2-body X + H problem. The H atom is supposed to be point-like 
and to interact with the incoming particle via a central potential 



V(r) 



1 a(r) 



(1) 



a(r) tends to the H dipole polarizability (a<2 = |) for large vàlues of r and 
regularizes the \ singularity at r=0. Its precise form is given in |J. 

We have displayed in Fig. [ï] the binding energies of the 7r + H bound states 
with interaction (jï|). One can see a large number of states with angular momen- 
tum vàlues up to L=7. Some of them are very close to the dissociation threshold. 
The elàstic /i + H cross section in the momentum range k £ [0.48, 0.64] has been 
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Figure 2. Elàstic /iH cross section for several 

Figure 1. Discrete spectrum for 7rH in the 

partial wave displaymg resonant benaviour 

two body approach 

plotted in Fig. Q with the most relevant partial wave contributions shown sep- 
arately. The L=5,6,7 states display a clear resonant behaviour also visible in 
the total cross section. A direct calculation of their position and width - us- 
ing the complex rotation method § - provides the vàlues E5=(7.9-1.4z)10~ 4 , 
E 6 =(l.l-0.25i)10- 3 , E 7 =1.0xl0" 3 -1.3ixl0" 7 respectively @. These resonances 
concern all high angular momentum states, involving large centrifugal barriers 
and relatively large energies. It's worth noticing however that they can be found 
even for L=l and at extremely small energy vàlues. In the 7r + H case one has 
e.g. Ei=(4.9-1.4i) x 10~ 7 . The examples shown in Figs. [ï] and ^ illustrate well 
the kind of physics governed by the polarization forces. An extensive work as a 
function of the projectile mass has been done || showing a very rich spectrum 
of bound and resonant states, with a complexity increasing as a function of the 
projectile mass. 

To what extend the results of this simple approach are reliable? Answer- 
ing this question was the main motivation of this work. A definite answer will 
only appear by letting the electron dynamics in H play its full role, that is by 
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considering the (X + p + e~) three body problem. And this is the aim of the next 
section. 



3 Three-body calculations 

The 3-body (X + p + e~) calculations are performed using Faddeev equations in 
configuration space. Three difïerent sets of Jacobi coordinates are involved, de- 
fined by 

1/2 

(r/j-r 7 ), (2) 



y< 



2m a {mp + m 7 ) 



1/2 



rpmp + m 7 r 7 

(»"« T — ) 

m/3 + m 7 
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where (a/?7) denote cyclic permutations of (123), and m a the particle masses. 
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Figure 3. Jacobi coordinates used in X + p + e~ calculations. 



We identify 1=X + , 2=p + , 3=e . The Standard Faddeev equations read 

(E-Hq- V a )V a = V a Y, Vfi, (4) 

where Ho is the 3-particle free hamiltonian and V a the 2-body Coulomb potential 
for the interacting (^7) pair 

V a (x a ) = . 6/367 . ej = ±l (5) 

For projectile masses mx < m p we restrict ourselves to scattering energies 
below the first rearrangement threshold X + +(p + e~)^p + +(X + e~). In that case, 
amplitudes $2 an d ^3 have no asymptotics. 

Equations (||) provide satisfactory solutions for bound states but are not 
suitable for scattering Coulomb problems. The reason is that their right hand side 
does not decrease fast enough to ensure the decoupling of Faddeev amplitudes 
in the asymptotic region and to allow unambiguous implementation of boundary 
conditions. In order to circumvent this problem, Merkuriev Q proposed to split 
the Coulomb potential V into two parts by means of some arbitrary cut-off 
function x 

V{x) = V s (x,y) + V l (x,y) 
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V s {x,y) = V(x)x{x,y) 
V l {x,y) = V(x)[l-x(x,y)] 

and to keep in the right hand side of equation (||) only the short range V s con- 
tribution. One is then let with a system of equivalent equations 

(E-H -W a - V£)9 a = V» Vp, (6) 

in which W a are some 3-body potential containing the long range parts: 

W a = VÏ + Vb + V* (7) 

This approach was found to be very efficient in calculating the e + Ps and e + H 
cross sections || [|. 




Figure 4. 3-body zero energy X H cross section as a function of the projectile mass mx 



Equations (|6|) were solved by expanding •f'j in the bipolar harmonies basis 



x%yi 



and their components <^j a4 in the basis of two-dimensional splines. In practical 
calculations we took the cutoff function: 





(aso) 




1 1 + exp 






1 + f 

yo 


j 



R. Lazauskas, J. Carbonell 



5 



Final results are independent of the parameters xo,yo,v but an appropriate 
choice for their vàlues makes the convergence of expansion (||) faster. The vàlues 
xq = 2.0, j/o ra 2Jmx, ^ = 2.3 are suitable. 

Binding energies for the lower 3-body 7r + H and /i + H bound states are plotted 
in Fig. ^. They are compared to the results of the 2-body approach (||). One 
can see that, although there is a qualitative agreement, 2-body energies are 
systematically underestimated. These results can be used to improve the short 
range part of the 2-body potential. 

Some interesting features of the 3-body Coulomb system can be learned from 
Fig. (Q), where the zero energy X + H cross section as a function of the projectile 
mass mx is displayed. Each peak corresponds to the appearence of a new S-wave 
bound state. The critical mass vàlues mi at which they occur, would enable 



to generalize the ground state stability triangle [ 10 1 to higher excitations. A 
zoom in the region of physical interest, /i + and 7r + , is shown in Fig. |5[ The 
calculated vàlues are respectively a^n = 69.1 and a n n = 24.4. Some care has to 
be taken in extracting the scattering observables specially at zero energy, from the 
asymptotic solution at finite distance. The long range polarization force makes 
the convergence of the observables as a function of the X + -H distance very slow 
and requires an appropriate extrapolation procedure M]. 
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Figure 5. Scattering length as a function of Figure 6. Energies for lower 7rH and /xH 
mx, in the region of fi and tt. bound states compared to 2-body results (|ï|) 

pH scattering deserves special comments. The 3-body wave function has to 
be antisymmetyric with respect to the p exchange. This can be realized in two 
different ways following the proton spin coupling. For the case when the two 
protons spins are antiparallel (singlet) the spatial part of the wave function is 
symmetric, while for the parallel case (triplet) it is antisymmetric. In the 2-body 
approach, these two cases give rise to completely different potentials. 

The singlet case has a broad attractive well which supports a great number 
of bound states. They have been calculated since the first days of Quantum 
Mechanics and they are presently known with a very high precision (see e.g. 



[11, 12 1 and reference therein). Our 3-body calculations cannot reach this kind 



of accuracy for bound states but are in good agreement for the lower excitations. 
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They provide furthermore the first result for the pH scattering length a s = —29.3. 
We notice that the zero energy scattering wave function shows 20 nodes in y\- 
direction, indicating the existence of 20 L=0 a g energy levels for ■ 

The triplet case - modeled by Landau 14] - is dominated by the Pauli 
repulsion between the two protons, overbalanced at r ~10 by the attractive 
polarization forces. 
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Figure 7. Zero energy pH in the pp triplet state, compared to the results of Landau potential 



Our 3-body calculations give a scattering length of at = 750 ± 5. The nodal 
structure of the Faddeev amplitudes indicates that such a big value is due to 
the existence of a first excited L=0 state with extremely small binding energy. 
By using the modified effective range theory [15| we are able to determine its 
binding energy, which turns to be B=(1.135±0.035)xl0 -9 below the H ground 
state. To our knowledge, this is the weakest bond ever predicted, three times 
smaller than the 4 He atòmic dimer [13|. 

The existence of L=0 and L=l ground states is well known and their binding 
energies have been very precisely calculated [|ÏÏ], 12]. These authors were however 
not able to conclude about the existence of a second S-wave bound state. 

It is interesting to compare the three-body calculations with those provided 
by the simple Landau two-body potential. S-wave cross sections are plotted in 
Fig. [7]. At zero energy, both calculations differ by two orders of magnitude while 
at energies E~10 -6 they are already in quite a good agreement, despite the 
simplicity of the 2-body approach. This is due to the fact that the effective 
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interaction is highly repulsive at short distances, what prevents the incoming 
proton to penetrate inside the H atom and minimizes the effect of the 3-body 
dynamics. 




E (u.a.) E (u.a.) 



Figure 8. pH elàstic cross section for L=3 in Figure 9. pH elàstic cross section for L=4 in 
the pp spin triplet state the pp spin triplet state 

The pH cross sections for higher partial waves have also been calculated. 
They exhibit some narrow resonances in several partial waves. Figs. | and | 
show the elàstic cross section for the L=3,4 states. They are compared to the 
results of Landau potential. We can see that their agreement in this energy 
region is rather good. The 3-body resonances are a little bit shifted to the lower 
energy region and have smaller widths. The same effect is seen in the bound 
state calculations where the 2-body results are always slightly underbound. The 
position and width of these resonances were estimated to E=(5.13-1.61í) x 1CP 6 
for L=3 and E=(1.56-0.94í) x 10~ 5 for L=4. 

4 Conclusions 

By solving the Faddeev equations in configuration space, we have calculated the 
bound and scattering states for the (X + p + e~) system with projectile masses 
mx larger than the electron. The long range polarization forces give rise to a 
rich spectrum of bound and resonant states with increasing complexity as mx 
increases. 

Predictions for the physical cases are obtained. We found 

in particular the scattering length vàlues a 7T + H =2AA and a fl + H =69.1 a.u. 

Of special interest is the pH system in the pp spin-triplet state. We predict 
a second S-wave bound state with binding energy B=1.14xl0~ 9 a.u. below the 
H ground state. This constitutes the weakest bond ever predicted, even smaller 
than the 4 He atòmic dimer. 

The existence of such a nearthreshold state dominates the low energy pH 
scattering and results into a scattering length of af=750±5 a.u. A low energy 
proton approaching an H atom will behave like colliding with a large nanoscopic 
object. Several resonances occurring in different partial waves, but visible in the 
total cross section, are also predicted. The experimental confirmation of these 
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results would be very interesting. 

The calculations presented here are performed using a fully non relativis- 
tic dynamics with Coulomb pair-wise interactions. In this framework, they are 
parameter free with the only input of partides masses and charges. The perturba- 
tions induced by strong X + -p interactions or higher electromagnètic corrections 
have not been included. 

In view of the extreme sensibility of the pH results, it is necessary to quantify 
both relativistic and strong-interaction effects. A direct measurement of the pH 
cross section at very low energy seems unlikely. One can however access the low 
energy pH continuum in the final state of the photodissociation cross section. 
Work is in progress in these two directions. 

Acknowledgement. The authors are sincerely grateful to C. Gignoux for useful discussions and 
helpfïïl advises. The numerical calculations were performed at CGCV (CEA Grenoble) and 
IDRIS (CNRS). We thank the staff members of these organizations for their constant support. 

References 

1. J. Carbonell, J., F. Ciesielski, F., Gignoux, C: Few-Body Systems Suppl. 8, 
428 (1995) 

2. Voronin, A., Carbonell, J.: Phys. Rev. A57 4335 (1998). 

3. Mott, N.F., Massey H.S.W.: The Theory of Atòmic Col·lisions. 
Oxford Science Publications 3 rc ^ Edition 1965 

4. Moiseyev, N.: Phys. Rep. 302 211 (1998) 

5. R. Lazauskas, R.: Rapport de Stage DEA, ISN-UJF (2000) 

6. R. Lazauskas, R., Carbonell, J.: to be published 

7. Merkuriev, S.P. : Ann. Phys. 130, 395 (1980) 

8. A. A. Kvitsinsky, J. Carbonell, C. Gignoux: Phys Rev. A46 (1992) 1310 

9. A. A. Kvitsinsky, J. Carbonell, C. Gignoux: Phys. Rev. A51 (1995) 2997 

10. A. Martin, J.-M. Richard and T.T. Wu, Phys. Rev. A46 (1992) 3697 

11. Taylor, J.M., Yan, Z., Dalgarno, A., Babb, J.F.: Molec. Phys. 97, 25 (1999) 

12. Hilico, L., Billy, N., Gremaud, B., Delande, D.: Eur. Phys. J. D12, 449 (2000) 

13. F. Luo, G. McBane, G. Kim, C. Giese, W. Gentry, J. Chem. Phys. 98 3564 
(1993) 

14. Landau, L., Lifshits, E.: Mecanique Quantique, Ed. Mir Moscou 1975, 

15. O'Malley T.F., Spruch L., Rosenberg L., J.: Phys. Rev. 125 491 (1961) 



